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This paper is one of a series of five papers in a special session organized by the NASA 
Fundamental Aeronautics Program that addresses uncertainty assessments for CFD sim- 
ulations in hypersonic flow. Simulations of a shock emanating from a compression corner 
and interacting with a fully developed turbulent boundary layer are evaluated herein. Mis- 
sion relevant conditions at Mach 7 and Mach 14 are defined for a pre-compression ramp 
of a scramjet powered vehicle. Three compression angles are defined — the smallest to 
avoid separation losses and the largest to force a separated flow engaging more compli- 
cated flow physics. The Baldwin-Lomax and the Cebeci-Smith algebraic models, the one- 
equation Spalart-Allmaras model with the Catrix-Aupoix compressibility modification and 
two-equation models including Menter SST, Wilcox k — ui 98, and Wilcox k — ui 06 turbu- 
lence models are evaluated. Each model is fully defined herein to preclude any ambiguity 
regarding model implementation. Comparisons are made to existing experimental data 
and Van Driest theory to provide preliminary assessment of model form uncertainty. A 
set of coarse grained uncertainty metrics are defined to capture essential differences among 
turbulence models. Except for the inability of algebraic models to converge for some sep- 
arated flows there is no clearly superior model as judged by these metrics. A preliminary 
metric for the numerical component of uncertainty in shock — turbulent-boundary-layer 
interactions at compression corners sufficiently steep to cause separation is defined as 55%. 
This value is a median of differences with experimental data averaged for peak pressure and 
heating and for extent of separation captured in new, grid-converged solutions presented 
here. This value is consistent with existing results in a literature review of hypersonic 
shock — turbulent-boundary-layer interactions by Roy and Blottner and with more recent 
computations of MacLean. 


Nomenclature 


Some variable names are specific to a particular turbulence model. If they are not found in this list they 
will be defined in the section they are encountered. 


Roman symbols 

A+ 

B q 

c 

Cp 

Cf 


modified Van Driest coefficient [Eq. 10] 
dimensionless heating, q w /(p w c p u T T w ) 
speed of sound [m/s] 

heat capacity at constant pressure [J/kg-K] 
skin friction coefficient, t w / (J^p e u 1 2 e ) 


* Senior Research Engineer, Aerothermodynamics Branch; AIAA Fellow 

* Senior Research Engineer, Aerothermodynamics Branch; AIAA Associate Fellow 
Senior Project Engineer, AIAA Member 


1 of 44 


American Institute of Aeronautics and Astronautics 



Roman symbols 

Cf 

D s 

D t 

e 

e* 

E 

T 

h 

H 

n 

Is 

k 

l 

M t 

M t 

M 


n 



P 

P+ 

V 

Pr t 


q 


Qrad 


'Tprodjim 

Reg 

Reg 

Set 



tij 

T 


u 

'U'i 5 'Uj 
U T 

u 

u+ 

w s 

x 

Xi , Xj 

y 

y + 


incompressible friction coefficient 
effective diffusion coefficient for species s [m 2 /s] 
turbulent diffusion coefficient [m 2 /s] 
static energy [J /kg] 

vibrational energy per unit mass evaluated at temperature T, [J /kg] 
metric of difference between computation and experiment 
vorticity function [Eq. 17] 
static enthalpy [J /kg] 
total enthalpy, [J /kg] 

Heaviside step function 

first ionization energy of species s [J /kg-mole] 

turbulent kinetic energy [J /kg] 

mixing length [Eq. 9] 

turbulence Mach number, \/2k/c 

Mach number based on friction velocity, u T /c w 

molecular weight [kg/kg-mole] 

coordinate orthogonal to wall [m] 

molar rate of production of species s by electron impact ionization [kg-mole/m-s] 
pu T n/ p, normalized distance to wall 
pressure [N/m 2 ] 

transformed pressure gradient, [Eq. 11] 

production term in turbulent kinetic energy equation [Eq. 27] 

turbulent Prandtl number 

heat transfer rate [W/m 2 ] 

radiative energy transfer rate [J /m 3 -s] 

maximum allowed value of production to destruction ratio [Eq. 27] 
momentum thickness Reynolds number 
incompressible momentum thickness Reynolds number 
turbulent Schmidt number, /x*/ pD t 
mean-strain-rate tensor, \ 

zero-trace version of mean-strain-rate tensor, Sij — 77 

c _ 1 du k x 
2 dx k 

molecular stress tensor, 2 /j,Sij 
temperature [K] 
velocity [m/s] 

velocity component in i and j directions, respectively [m/s] 

friction velocity, s/tw Ip w 

velocity tangent to wall, [m/s] 

u/u T1 dimensionless velocity 

mass rate of production of species s [kg/m 3 -s] 

distance along wall (flat plate), coordinate in streamwise direction [m] 
coordinates in i and j directions, respectively [m] 
distance normal to wall (flat plate), coordinate orthogonal to x [m] 
pu T y/n , normalized distance (flat plate) 
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Greek symbols 


7 i 

intermittency function, Eq. 15] 

5 

boundary layer thickness, [m] 

5k 

velocity thickness, [Eq. 14] 

e+ 

eddy viscosity, pt/p 

4 

inner layer eddy viscosity 

4 

outer layer eddy viscosity 

7? 

molecular translational-rotational conductivity [J /K-m-s 

Vt 

turbulent translational-rotational conductivity [J /K-m-s] 

Vv 

molecular vibrational-electronic conductivity [J /K-m-s] 

Pvt 

turbulent vibrational-electronic conductivity [J/K-m-s] 

p 

viscosity [kg/m-s] 

V 

eddy viscosity, p/p 

P 

density [kg/m 3 ] 

T 

shear [N/m 2 ] 

T~ij 

Reynolds stress tensor [Eq. 28] 

Xs 

mole fraction of species s 


f duk-i du k+ 1 \ 

^3a;fc + i dxk-i ) 

Subscripts 


e 

at edge of boundary layer 

e— 

electron 

i 

component in i direction 

j 

component in j direction 

mol. 

molecules 

s 

species index 

t 

turbulent value 

V 

vibrational-electronic energy component 

VD 

parameter used in Van Driest transform 

w 

conditions at wall 

00 

reference condition in free stream 


I. Introduction 

The Fundamental Aero Program (FAP) has a goal of reducing uncertainties in hypersonic flow simulation 
by 50%. The target goal started as a subjective assessment of the simulation needs by NASA researchers 
looking at diverse challenges and rolled up into a single metric. The FAP has now been tasked to define the 
metrics by which success or failure is judged. 

The definition of simulation uncertainty in hypersonics has little meaning without specifying free stream 
conditions and configuration geometry. Different physical phenomena engage as a function of total enthalpy, 
atmospheric composition, vehicle size, and vehicle shape. Consequently, five mission relevant simulations 
have been defined to characterize simulation uncertainties for specific conditions. The focus here is on 
hypersonic cruise conditions for a scram-jet with shock wave - turbulent boundary layer interaction at a 
compression corner. We include compression corner angles that are sufficiently steep to cause separation of 
an incoming turbulent boundary layer (something not wanted for a scramjet) because it provides significant 
challenges to computational simulation. Furthermore, the ability to predict the maximum compression angle 
that does not separate the boundary layer provides opportunity to open design space and shorten the running 
length to the inlet. 

Quantifying uncertainty reduction further requires a well-defined initial simulation capability from which 
new modeling improvements can be measured. Consequently, we focus on turbulence models that have seen 
wide application in production codes through 2006 — the approximate time frame when FAP goals were first 
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announced. By “production” codes we refer to codes like LAURA, 1 ’ 2 DPLR, 3 and Vulcan 4 that are routinely 
used for simulation of hypersonic flows within NASA. The models used in these codes in the 2006 time frame 
include algebraic, and various one- and two-equation turbulence models. Detached eddy simulations (DES), 5 
large eddy simulations (LES), 6 and direct numerical simulations (DNS)' in the hypersonic domain are not 
considered in the 2006 baseline. Advances were clearly being made in these approaches 8 in that time frame 
but they were not yet routine or robust components of production codes for hypersonic simulations. In like 
manner, these production codes did not have routine capability to predict transition location using a coupled 
stability analysis like STABL 9 or LASTRAC.3d. 10 Specification of transition location, if implemented at 
all, generally was based on R eg/M e correlations. The dependence of the shock wave turbulent boundary 
layer interaction (SWTBLI) at a compression corner as a function of distance from the upstream transition 
location is not well understood. Presumably, if the transition location (tripped or natural) is close to the 
SWTBLI then the unsteady nature of transition will feed an additional source of unsteadiness into the 
interaction. A transition location far upstream (hundreds of boundary-layer thicknesses) feeding a fully 
developed, turbulent boundary layer profile into the SWTBLI is therefore a prerequisite for experimental 
data sets used to establish uncertainty. Even with a fully developed, incoming turbulent boundary layer 
experimental results show that shock oscillation is of inherent nature in the SWTBLI with separation and 
that shock oscillation is considered to be the consequence of the coherent structures in the separated region. 11 

This uncertainty definition process is started by building on the 2006 review of turbulence models for 
hypersonic flows by Roy and Blottner, 12 which in turn sought to update the comprehensive review of hyper- 
sonic shock/turbulent boundary-layer interaction experiments published in 1991 by Settles and Dodson. 13,14 
Knight et al. 15 have assembled a complementary review for supersonic flow including interactions associated 
with compression corners, shock impingement, single fin, and double fins.; 15 however, the Mach numbers 
(2.3, 3) and adiabatic boundary conditions for the compression corner problems do not address the hyper- 
sonic domain, especially the role of a cold wall boundary. More recent experiments with computation by 
Reinartz et al. 16 examined wall temperature effect in a shock tunnel flow at Mach 7.7 and 7.4 but did not 
have fully developed turbulent flow on the flat plate leading to the 15 degree compression corner. Holden 
et.al. 1 ' have provided new contributions to the computational and experimental analysis of the canonical 
SWTBLI compression corner problems in the LENS facility at CUBRC. These new data sets are used herein 
because they feature long models with natural transition far upstream of the compression corner at hyper- 
sonic conditions and the compression surface is sufficiently long to allow a plateau in pressure downstream 
of the separation without a terminating expansion to suppress growth of the recirculation. 

The paper is structured as follows. We first define a mission relevant test problem including free stream 
conditions for which the simulation uncertainty will be characterized. The governing equations and tur- 
bulence models are then documented to remove any ambiguity of how they are formulated. A suite of 
uncertainty quantification metrics are introduced based on a collection of simulations reviewed by Roy and 
Blottner. 12 This suite of coarse grained uncertainty metrics is defined to easily capture their values from 
figures in the literature. Relevant experimental data on compression corners are then reviewed and new 
simulations are executed with a variety of models and codes representative of the state of the art in 2006 as 
explained earlier. The new metrics are calculated for both existing and new simulations and a representa- 
tive value is proposed. Finally, we return to the mission relevant test problem (which has no experimental 
reference data), execute new simulations including gas chemistry perturbations in order to assess an overall 
uncertainty metric. Two appendices are also provided to present verification and validation data on a simpler 
flat plate problem. 


II. Definition of Mission Relevant Test Problem 

The development of a hypersonic cruise vehicle using air breathing propulsion is deemed a high pri- 
ority mission within NASA’s Fundamental Aero Program (FAP). Two trajectory points are defined for a 
representative mission 18 in Table 1 below. The geometry for the compression ramp leading to the inlet is 
a two-dimensional, double-wedge approximation (Fig. 1) that is loosely based on an X-43A configuration 
discussed by Berry et al. 19 The two-dimensional configuration has a 4.5 degree wedge with a 12 ft. (3.66 m) 
running length from the sharp leading edge to the compression corner. The second section provides another 
5.5 deg. compression (10 deg. total) with another 6 feet (1.83 m) running length to the end of the body. (A 
blunted leading edge with 1 cm radius was rejected in favor of the sharp leading edge; it exhibited entropy 
effects that influenced laminar separation at the corner.) 
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Table 1. Trajectory Points for Mission Relevant Problem 



Too, m/s 

Poo, kg/m 3 

alt., km 

T, K 

T W all, K 

7 

2086 

4.401 10“ 2 

24.4 

220.94 

500 

14 

4418 

9.811 10“ 3 

33.7 

246.93 

500 



Figure 1. Pressure contours over mission relevant configuration, (blue low, red high) 


The 5.5 deg. compression is not intended to separate the flow at the corner. Two additional configurations 
with a 30 deg. compression (34.5 deg. total) and a 35 deg. compression (39.5 deg. total) will also be analyzed. 
In these cases, a separation zone is desired because simulation of a separated zone is more challenging and 
because the ability to identify an incipient separation condition is important for understanding design space 
constraints. These configurations and conditions are deemed canonical problems with relevance to a high 
priority mission of interest to FAP. 


III. Conservation Equations 

The conservation equations with effects of RANS turbulent dissipation are provided below. Turbulent 
kinetic energy, if modeled by a 2-eq. formulation, is included in the total energy conservation equation. Note 
that source terms for chemical reaction and thermal relaxation do not include any modification associated 
with real turbulent fluctuations. 


A. Species Conservation 

d(p a ) , d(p s Uj) _ d 


dt 


dx-j 


dxi 


pD s 


p t M s \ d\ 


S Ct M J dxj j 


■ w s 


(1) 


B. Momentum Conservation 

d(pui) d(pujUi) dp d 


dt 


+ 


dXn 


~ dx t + dx 5 [tji + Tji] 


(2) 


C. Total Energy Conservation 


d r / , UiUi 

p[e+ — + k 


dt 


where 


d 

dxj 

d 

dxj 

d 

dxi 


puj (h + ^ + fc) 


. , .. d 

[Ui ( tji + Tji)\ + — 




( 3 ) 


, dT dT v 

{ri + Vij-K— + {Vv + Vvt)-K— 

(J JL j (J Jb j 


d 

dxi 


pk\ dk 

P + CTk 

u> / dx-i 


Qrad 


V th Vv Pt 

It = — ; t—, p vt = 


P + Vv c p Pr t 


V + Pv c p Pr t 


( 4 ) 
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D. Vibrational-Electronic Energy Conservation 


dpe v 

dt 


dpu^e v du.j 


d 

dx-j 


E ? D 


j 

S^M 



d ' 

1 

B 

+ dxj 


, ^ sdTv 

{Vv + Vvt)-Q^ 


J2 Ps (e ^ e »■«) + 2 Pe _ | R(T - T v ) Y, n e -J. + E ™sD s -Q rad 


s—mol. 


s^e— 


s—mol. 


( 5 ) 


(6) 


IV. Turbulence Models 

The turbulence models used herein are fully documented in this section to avoid any ambiguity in the 
value of constants or the definition of source terms. In this regard, a current, on-line resource on turbulence 
models may be found at http://turbmodels.larc.nasa.gov. 


A. Algebraic Models 

The algebraic formulations of eddy viscosity use the Gupta modifications 20 to the models of Cebeci-Smith 21 
and Baldwin-Lomax. 22 The model formulations are repeated below to maintain a complete record of all 
turbulence models used within LAURA 1 for this study. The original references provide a more complete 
explanation of the model derivations. 

The algebraic eddy viscosity is given by 


e 


+ 


77. S 7 ^crossover 
77 77 crossover 


where n crossover is the value of n where ef = e+. 


( 7 ) 


1. Cebeci-Smith 


The inner eddy viscosity in the baseline model is defined as 


e+ 


A + 


P+ 


pl 2 dU 
fi dn 

0.4rr [l — exp (— n + /A + )] 
26 

vT - 11.8P+ 

_p_ / dp\ 

p 2 v? T \<9s/ e 


(8) 

(9) 

( 10 ) 

( 11 ) 


An alternative definition of A + suggested by Gupta 20 and found to yield good agreement with turbulent 
heating data on blunt bodies for cold hypersonic flow is 


A + 


26 


(12) 


The approximate equivalence of t/t w and 1 — 11.8P + can be derived by numerically integrating the stream- 
wise momentum equation across the sublayer from y + = 0 to y + = 11.8 as discussed by Cebeci. 23 The 
definition of A + in Eq. 10 removes the singularity that can occur in the case of separated flow not to men- 
tion the implicit dependence of r on ef . However, strongly expanding flows (as around the corner of a blunt 
body) can lead to negative values of the argument in Eq. 10. Results presented herein use Eq. 12 unless 
otherwise noted. 
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The outer eddy viscosity is 


e 


+ 

O 


Sk 


7 i 


0. 0168 pUeSkji 



1 + 5.5 



(13) 

(14) 

(15) 


The coordinates n and s are orthogonal and tangent to the wall and U is the velocity component in the s 
direction. The boundary layer thickness 5 is defined as the value of n where = 0.995. This definition 
needs to be adjusted in the case of shock layer radiation where the total enthalpy is not conserved across 
the inviscicl layer. 


2. Baldwin-Lomax 


The inner eddy viscosity is computed as in the Cebeci-Smith model using Equations 8, 9, and 12. The outer 
eddy viscosity is 


_ 


Hn) = 


0.0l68CcPpF wa keFKleb 


P 


n 


dU 


dn 


[l — exp (— n + /A + )~\ 


F> 


Kleb — 


1 + 5.5 


0.3?i 


-l 


(16) 

(17) 

(18) 

(19) 


Here Ccp = 1.6, T max is the maximum value of the vorticity function T across the outer layer and n max is 
the value of n where F(n max ) = T max ■ Note that the thin-layer approximation |^| ss |fi| is used to define 


B. One-Equation Model — Spalart-Allmaras 

The compressible form of the Spalart-Allmaras 24 model developed by Catrix and Aupoix 25 is defined 


d(pP) d(puji>) 


dt 


dx 


Id ( dv \ 1 d ( d^fpv \ c b2 d^/pp dy/p, 


J ' = CbiHpD-c^i/^p ( - ) +--^~ ( P^ ) +-tt- ( \Tpv~ ~ 

' a J a axd \ oxd J a oxn \ axn 


a dxj dxi 


(20) 


The turbulent eddy viscosity is defined 


v 

Pt = fvl = — 5 . g - ) Xv — 

Xt + c 4 i v 

Other constants and functions required to complete this model are: 


(21) 


q,i = 0.1355, o >2 = 0.622, c„i = 7.1, a = -, k = 0.41 

o 

_ 0>1 1 + 0)2 no „ _ o 

CujI — 9 H - C-cu2 — U.o, Ccj3 — 2 


(22) 

(23) 


U =9 


1 + c w3 

1/6 

g = r = c^ii" 6 — r), r = min 

-^—,10 

_9 6 + c ®3- 



[Sn 2 d 2 J 


’ ( 24 ) 
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C. Two-Equation Models 


The two-equation models used herein share the following conservation law form for turbulent kinetic energy 
k and the specific dissipation rate w . 26 


d(pk) | d(piijk) 
dt ■ U 
d{pu) 
dt 


dccj 


= 'P ~ (3*puk + 

d 


P 


P 


d(puju) 7 

C j 1/ KJ du j 

All of the two-equation models considered herein define the production term V as 


= -P- PcPCO 2 + 

dxn v* dx 


pk\ 

dk 


dxj. 

pk\ 

duj 


dxj. 


■ i dui * 

V = mm ( Tij , r prod ,iimPc P^k 

2 

2//f S 7 j —pkdij 


, dk du 
dxj dxj 


(25) 

(26) 

(27) 

(28) 


The production term is limited by some factor r pro( ]j irn (defaults to 20 in LAURA) times the destruction 
term to provide more realistic values of turbulent kinetic energy behind shocks. 


1 . Compressibility Corrections 

Rumsey 27 has reviewed the impact of compressibility corrections on turbulent, hypersonic boundary layers 
on flat plates. He investigates the claim that “... compressibility corrections are not required for hypersonic 
boundary layer flows ... for a wide range of Mach numbers and wall-temperature boundary conditions” and 
ultimately notes that compressibility corrections derived for free shear layers do not work well for boundary 
layers and vice versa. Compressibility corrections were originally derived to account for the reduction in 
a free shear layer growth due to compressibility . 28 In dealing with separated flows, it seems logical that 
some type of compressibility correction is needed though an optimal implementation is not obvious. The 
Sarkar/Zeman 29,30 compressibility modifications to define /3* and (3 C from the baseline model values of (3* 
and (3 are applied here using Wilcox’s model to define F(M t ). 

P* c =P*[l + tf(M t )}, 0 c = 0-0TF{M t ) (29) 

with 

r = 2, M t0 = i, F(M t ) = [M 2 — M 2 o ] H(M t — M to ) (30) 

In general, compressibility corrections are not applied (. F(M t ) = 0) unless specifically noted. 


2. k-w 2006 26 

The turbulent eddy viscosity is defined as 

pk 


Pt = 


u = max 


ui, Cii . 


'2 SijSij 


(3* 


Clim — 0 


Other constants and functions required to complete this model are: 

13 

(Tfc = 0.6, cr w = 0.5, 7 = — 


U = 


1 + 85 X „ 


Xu = 


1 + 100y ( 

Ghost cell values for k and ui behind a no-slip wall are defined as 

z - n _ 617 

k gc - 0 , u gc - pod2 


(31) 


* 

V - 

k 

= -, /?* = 0.09, /3q = 0.0708, 

UJ 

P = Pofp, 

(32) 

kj $ki 

1 _ P a d _ j 0 for 

7 0 1 1 1 r 

w y ^ for 

dk du) rv 

dxj dxj — 

(33) 

to ) 3 

dk du ^ rj 
dxj dxj 


(34) 


where subscript 1 indicates the value at the bounding cell center and d is the distance from the wall to the 
bounding cell center. 
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3. k-u 1998 26 


The turbulent eddy viscosity is defined as 


Mt = 


pk 

U! 


Other constants and functions required to complete this model are: 

07 = 0 . 5 , a u = 0 . 5 , 7 =^> 0 = 0 

25 id 


A)= 0.0708, /3 = 0ofp, ft i = 


1 


00 — 0.09, j3* — Pof/3*, f/3* — < l+680xfc 

l l+400 X fc 





= 

Xk 

< 

0 




0 ’ 

Xk = 

Xk 

> 


as 





&ij ^kj &ki 

(ft*) 3 

1 dk dm 
lu 3 dxj dxj 


, __ n 6/q 

fc gc - 0, u; gc - 


(35) 


(36) 


(37) 


where subscript 1 indicates the value at the bounding cell center and d is the distance from the wall to the 
bounding cell center. 


4- Menter-SST 31 

The turbulent eddy viscosity is defined as 


Mt = 


pa\k 


a x = 0.31 


max(aiw, O F 2 ) 

Other constants and functions required to complete this model are: 

\fk 500^ \ Apa u2 k 


(38) 


Fi = tanh.(arg\), argi = min 


max 


P*wd’ d 2 io ] ’ CD ku d 2 
F 2 = tanh(argf), argi = max (2 


, CD kuJ = max 

Vk 500 ^\ 


/ 2pa u2 dk dio 
\ to dxj dxj 


, 10“ 20 


p*u)d’ d 2 u 


(39) 

(40) 


<r k = F ± a kl + (1 - Fx)a k2 , a kl = 0.85, a k2 = 1.0 

07 = FiVui + (1 - Fi)a u2 , 07i = 0.5, <r u2 = 0.856 

(3 = Fxl \ + (1 - Fl)/3 2 , Px = 0.075, p 2 = 0.0828 

7 = F171 + (1 - Fi)72, 71 = Px/P* - Vu,xK 2 /\fp*, 72 = Pit P* ~ 072K 2 / \/P* 

P* = 0.09, v* = v t , 0 = 2(1 -Fx)^Fl, k = 0.41 

id 

Ghost cell values for k and u> behind a no-slip wall are defined as 

60iq 


kgc 0 , LO gc 


Pid 2 


(41) 


(42) 


where subscript 1 indicates the value at the bounding cell center and d is the distance from the wall to the 
bounding cell center. 
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V. Uncertainty Metrics / System Response Quantities 


A. Simulation Uncertainty 

As will be seen in the next section, a significant variation in computed results is observed as a function of 
turbulence model for shock wave / turbulent boundary layer interactions (SWTBLI) that induce separa- 
tion. Simple, coarse grain uncertainty metrics for system response quantities (SRQ in the terminology of 
Oberkampf and Roy 32 ) are proposed here to capture essential differences between computation and experi- 
ment. As additional experimental data is gathered and turbulence models are winnowed and converge, finer 
grain metrics are required to finish the evolution of turbulence models for hypersonic flows. 

The coarse grain metrics focus on surface measurements in the the pre-compression region, separated 
region (if it exists experimentally), and the post-compression region including the maximum measured values 
of surface pressure, heating, and shear. Off-body measurements are deferred for use in defining the fine grain 
metrics. The metrics in this section are computed without any regard for experimental uncertainty. The 
increments associated with experimental uncertainty will be computed in the next sub-section. 

In the post-compression region (see for example the region where x > 1.04 in Fig. 8) metrics are proposed 
that compare the maximum measured values of pressure, heating and shear to the maximum computed values. 
The maximum values need not occur at the same physical location. Thus, 


t-,/ r \ I Jpost.cfd.max J post. exp, max \ 1 AA 

E(fmax) = ~ ^“7 — X 100 

\ Jpost,exp,max / 


(43) 


where / represents p, q, or r. If there is no local maximum in the post-compression region the values 
immediately following the jump will be used. 

The metric in the pre-compression region (see for example the region where x < 0.95 in Fig. 8) is tuned 
to check that inflow conditions setting up the interaction are correctly simulated. Select a subset of N pre 
measurements just upstream of the interaction to form an average value of f pr e,exp and an average value of 
location Xpyg^gxp 


%pre,exp 


N v 


N„ 


^ f exp (' ! 'ri,exp ) 


N„ 


N pre 

^ ) Xn,exp 
n—1 


(44) 

(45) 


The subjective nature of choosing N pre upstream measurements is acknowledged. The definition is in keeping 
with the coarse grain metric definitions that are intended to be easily defined by inspection of plotted 
results. If there is a significant difference between the measurements that make up the average relative to 
the measurement uncertainty in this domain then this metric should be abandoned. Assuming an appropriate 
average value can be defined then the pre-compression metrics are given by 

E(fpre) = ^ W/«i-/ea: P (a: P re,exp) \ x 1QQ ( 46 ) 

\ Jpre,exp / 

The metric in the separation region (see for example the region where 0.95 < x < 1.04 in Fig. 8) is 
defined at the midpoint between the experimentally measured separation point and the corner location, 

• Vsep,mid (*E sep,exp T ^corner) /2. 


E(fsep) 


fcfd{x sep,mid) f exp(x sep,mid) 

f exp^X sep,mid) 


x 100 


(47) 


Note that the separation metric is defined even if the simulation shows no separation. Finally, a metric for 
the extent of separation is defined 


E(L) = f Xs£ P^ d X ^P.exp \ x 10Q ( 48 ) 

\ *Esep,exp % corner / 

The separation metric is defined at the location where the pressure and heating show an abrupt increase 
in slope. All of these metrics return negative values if the computed variable under-predicts the measured 
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variable. The value of E(L) is 100 if the experiment shows any separation and the computation shows 
no separation. Many of the metrics recorded here were extracted from published figures where there was 
difficulty in reading differences to better than a couple of percent. For this reason, all metrics are consistently 
rounded to the nearest 5%. 

A single metric may be defined from these coarse grain definitions to provide a simple measure of the 
current state-of-the-art. This metric is strictly used to provide an easily defined target (though still a 
challenging goal) for the Fundamental Aeronautics Program (FAP) to assess improvements in hypersonic 
simulation capability. This metric will be introduced after reviewing the results of multiple simulations here 
and in the literature. 

B. Experimental Uncertainty 

The CUBRC experimental data 17,33 is the prime source of experimental data here because the large models 
and test conditions produce natural transition well upstream of the interaction region. Of all the available 
experimental data sets, the CUBRC data best covers the Mach number range and total enthalpy range of the 
mission relevant problem as shown in Fig. 2. Uncertainty in measured pressure is given as ±3%. Uncertainty 
in measured heating is given as ±5%. Uncertainty in measured shear was not discussed. Other sources 34,35 
estimate shear measurement uncertainty at ±7% to ±12% with ±10% assumed for subsequent calculations. 
Uncertainties in free stream dynamic pressure (p 00 V 2 00 /2) and stagnation point enthalpy (sa V 2 oc /2) are 
±5%. Using dimensional analysis, the uncertainty in free stream energy flux (pooU 3 oo/2) is 1.05 x 1.05 0 5 — 1 
or ±7.6%. If one assumes that the changes in pressure coefficient, skin friction coefficient, and heat transfer 
coefficient at any point on the body are small compared to changes in these free stream quantities then 
the total experimental uncertainties due to calibration and free stream uncertainties are estimated to be: 
A p exp = V3 2 ± 5 2 = ±5.8%, A r exp = V10 2 ± 5 2 = ±11.2%, and A q exp = V5 2 ± 7.6 2 = ±9.1%. The 
uncertainty for separation length A L exp is taken as the distance between gauges at the separation point 
divided by the experimental separation length times 100. For Run 54 this uncertainty is computed as 
A L exp = (0.32/2.54)100 = 12.6%. These estimates do not include measurement errors associated with high 
frequency unsteadiness relative to gauge response time — possibly an issue inside the separation zone. Nor 
do they include any errors associated with flow non-uniformity in time and space or model alignment. 

E+06 E+06 




(a) Reynolds number versus Mach number (b) Reynolds number versus total enthalpy 


Figure 2. Relation of experimental data sets (symbols) to mission relevant problem trajectory points (con- 
nected by line). 


Some other experimental cases covered in the review by Roy and Blottner are also presented to provide 
a more diverse set of simulations from which uncertainty metrics may be derived. We assume experimental 
errors in these sources are at least as large as those discussed above from CUBRC. A more complete discussion 
of the definition of uncertainty metrics follows in the next section, after having the opportunity to observe 
the primary ways in which the simulation and experiment differ. 
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VI. Test Cases with Experimental Data for Compression Corners 


In the following figures the flow is in the x — y plane defined by z = 0. The x coordinate is aligned with 
the flow direction in the free stream. The wall shear stress is computed as r = cos (9)t x + sm(6)T y where 6 
is the inclination angle of the wall relative to the x axis. 

In all cases the turbulence model is engaged from the leading edge of the model even if the transition 
location is defined in the experimental data. There is an implicit assumption that fully developed turbulent 
flow is present upstream of the interaction region in all of the simulations. The possible effects of transition 
location and transition length are accepted as contributions to the modeling uncertainty. 

A. CUBRC Compression Corner - Run 54 

Experimental data for a turbulent boundary layer at Mach 11.3 on a flat plate followed by a 36 deg. ramp 
obtained at CUBRC 17 are used as the first validation check for the suite of turbulence models applied to 
shock / turbulent boundary layer interactions at compression corners. The geometry is presented in Fig. 3. 
The baseline grid uses 264 cells along the plate and 94 cells normal to the plate. Two additional grids are 
also utilized which provide a factor of two (fine grid) and 4 (very fine grid) additional resolution in both 
coordinate directions relative to the baseline grid. The near wall region prior to separation is resolved with 
y + « .015, achieved by specifying a cell Reynolds number equal to 0.1 in the initial grid generation process. 
Inflow boundary conditions are: = l,769m/s, = 0.08246kg/m 3 , T\ ' x = 61K and a = 0 cleg.. These 

inflow conditions are very close to those of the flat plate case used as a verification and validation check 
in Appendix B. Surface boundary conditions are no-slip with T w = 300K. Supersonic outflow boundary 
conditions are extrapolated from the interior. The test gas, air, is treated as calorically perfect unless 
otherwise noted. Molecular transport properties are computed from Sutherlands law. 



Figure 3. Turbulent wedge geometry for CUBRC Run 54. 17 


Figure 4 compares experimental data for pressure and heating along the surface to simulations using 
various turbulence models with LAURA. The solid line result in each figure shows the simulation and 
symbols show the experiment. A log scale is used to provide equivalent resolution of pre- and post- shock 
levels of heating and pressure. In general the percent difference between simulated and measured heating in 
this compression corner case is equivalent to that same difference measured for the flat plate in Appendix 
B although the log scale makes comparisons look better. For example, in the baseline k — to (2006) in 
Fig. 4 (d), the average measured heating on gauges with 0.8 < x < 0.96 is 4.86 W/cm 2 . The average 
simulation value over this same range is 6.43 W/cm 2 , an over-prediction of 32%. For the flat plate the 
average measured heating over the same range is 5.40 W/cm 2 and the average simulation value is 7.18 
W/cm 2 , an over-prediction of 33%. The coarse grained metrics from Eqs. 43-48 for this figure are captured 
in Table 2. 

All of the algebraic models failed to converge for this case. The separation point moved far upstream 
and large scale, unsteady motion ensued within the separation zone. The Cebeci-Smith model using the 
damping coefficient given by either Eq. 10 or Eq. 12 behaved similarly. The experimental data suggests ex- 
istence of a rather flat heating plateau. The remaining 2-eq. models under-predict the extent of separation 
if no compressibility correction is applied. A general trend is observed for the simulations without com- 
pressibility correction (Fig. 4 (b-d)) such that larger predicted separation zones produce greater overshoots 
in post-shock peak heating rate and, to a lesser extent, pressure. Heating overshoots in these cases vary 
from approximately 100% to 25%. The downstream asymptotic heating levels never recover with difference 
varying from approximately 100% to 20%. Application of a compressibility correction (with M t o > 0.4 to 
avoid conditions that disturb agreement with the van Driest distributions noted in Appendix B) produce 
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a separation extent in better agreement with experiment — see Fig. 4 (e-f). Pressure levels across the 
separation zone are in good agreement with measured data. Predicted heating levels in the separation zone 
are lower than measurements and no heating plateau is observed. 

Comparisons to experimental data for shear are presented in Fig. 5 using a linear scale to capture positive 
and negative values. In addition to over-predicting extent of separation as noted earlier the Spalart-Allmaras 
model under-predicts peak shear by 150%. The Menter SST model has best agreement with the shear data 
on the baseline grid. All of the models tend to under-predict the measured shear while over-predicting the 
measured heating in the post-compression domain. 




(a) l-eq: Spalart-Allmaras (b) 2-eq.: Menter SST 




(c) 2-eq.: k — u (1998) (d) 2-eq: k — u> (2006) 




(e) 2-eq: k - uj (2006) (CC M t0 = 0.5) 


(f) 2-eq: k - uj (2006) (CC M t0 = 0.4) 


Figure 4. Pressure and heating distributions over a compression corner at Mach 11.3 computed with 1- and 
2-eq. turbulence models on baseline grid (264 x 94). Algebraic models failed to converge for this case. 


All of these compression corner results were computed on grids that yield grid converged results for the 
flat plate. However, the resolution of a separated region prior to a compression corner generally requires 
finer resolution. A previous validation study for laminar flows 36 shows that onset of separation tends to be 
under-predicted (smaller recirculation) until a grid converged solution is achieved. Based on this observation, 
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the k — u (2006) model with various compressibility corrections is retested on finer grids (fine: (528 x 188) 
and very fine (1056 x 376)). The recirculation region is resolved with approximately (148 x 228) cells on the 
very fine grid. 

The surface pressure and heating distributions on these uniformly refined grids are presented in Fig. 6 and 
the surface shear distributions are presented in Fig. 7. The simulations without compressibility correction 
fail to show any significant separation even on the finest grid. The simulation using M t o = 0.7 provides the 
best overall agreement with measured pressure, heating, and shear levels (Fig. 6 (b,f) and Fig. 7 (b,f)). 
The extent of separation appears to be grid converged. The heating rate overshoot on the post compression 
side is approximately 130%. Some additional structure appears in the heating distribution on the finest grid 
that was not evident on the fine grid. The fine grid heating in the plateau region is in good agreement with 
measurement but an approximately 70% overshoot is evident on the finest grid. The coarse grained metrics 
from Eqs. 43 - 48 for Fig. 6 are captured in Table 3. 

Table 2. Coarse-Grained Uncertainty Metrics for Fig. 4 


Model 

E {c[max ) 

E(jpmax) 

E(qsep) 

E(p S ep') 

E^Qpre') 

E(jppre) 

E(L) 

SA, (264 x 94) 

0 

40 

50 

20 

40 

10 

20 

SST, (264 x 94) 

100 

35 

35 

10 

15 

10 

20 

k-uj 98, (264 x 94) 

60 

15 

-220 

-no 

15 

10 

80 

k-w 06, (264 x 94) 

20 

10 

-230 

-90 

20 

10 

95 

k — u) 06.5, (264 x 94) 

70 

40 

-100 

5 

20 

10 

0 

k — ui 06.4, (264 x 94) 

60 

50 

-100 

-25 

10 

10 

-5 


Table 3. Coarse-Grained Uncertainty Metrics for Fig. 6 and associated, un-plotted results from Vulcan and 
DPLR 


Code 

Model 

E(qmax ) 

E{jPmax ) 

E(q S ep) 

E(p S ep ) 

E(qpre') 

E(j)pre) 

E(L) 

LAURA 

k — u) 06, (528 x 188) 

20 

10 

-230 

-90 

15 

10 

100 

LAURA 

k — ui 06.7, (528 x 188) 

120 

80 

0 

40 

15 

10 

0 

LAURA 

k-u 06.6, (528 x 188) 

110 

80 

-30 

40 

15 

10 

-10 

LAURA 

k — u 06.5, (528 x 188) 

110 

90 

-60 

-30 

15 

10 

-60 

LAURA 

k-uj 06, (1056x376) 

25 

10 

-230 

-90 

15 

10 

100 

LAURA 

k-u 06.7, (1056x376) 

130 

80 

30 

50 

15 

10 

0 

Vulcan 

k-uj 06, (1056x376) 

60 

35 

30 

35 

15 

10 

-50 

Vulcan 

k — u 06, uni 

60 

40 

40 

10 

15 

10 

0 

Vulcan 

k-u 98, (1056x376) 

60 

50 

40 

55 

15 

10 

-55 

Vulcan 

k — ui 98, uni 

5 

15 

-70 

-155 

15 

10 

100 

Vulcan 

Men. SST, (1056x376) 

70 

60 

60 

70 

15 

10 

-190 

Vulcan 

Men. SST, uni 

75 

60 

60 

60 

15 

10 

-140 

DPLR 

Men. SST, (266 x 96) 

40 

45 

10 

0 

15 

10 

-35 

DPLR 

Men. SST, mod. stress lim. 

45 

35 

0 

-5 

35 

10 

5 


In Fig. 8 we focus on a single turbulence model (Menter - SST) but consider how the simulation results 
may vary as a result of different implementations within a code (Fig. 8-a) or between codes (Fig. 8-b). The 
progression of models study was motivated by an observation that the extent of separation observed in DPLR 
results (V4.01.1) for the Menter-SST was over-predicted compared to experimental data 17 but was under- 
predicted in the LAURA implementation (Fig. 4-b). (See a related study in Appendix B for a flat plate.) 
The DPLR solution was rerun using the LAURA grid and perfect-gas model. 37 No significant change in the 
extent of separation was observed. The remaining steps made changes to the LAURA implementation to 
make it more like the DPLR implementation. The DPLR implementation of the conservation laws for k and 
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(a) 1-eq: Spalart-Allmaras 



(b) 2-eq.: Menter SST 





(c) 2-eq.: k — co (1998) 





(d) 2-eq: k — to (2006) 




(e) 2-eq: k - co (2006) (CC M t0 = 0.5) (f) 2-eq: k — lj (2006) (CC M t0 = 0.4) 

Figure 5. Shear force distributions over a compression corner at Mach 11.3 computed with 1- and 2-eq. 
turbulence models on baseline grid (264 x 94). Algebraic models failed to converge for this case while showing 
massive, unsteady separation. 
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(a) No compressibility correction, on fine grid (528 (b) CC Mto = 0.7, on fine grid (528 x 188) 

x 188) 




(c) CC Mto — 0-6, on fine grid (528 x 188) (d) CC Mto — 0-5, on fine grid (528 x 188) 




(e) No compressibility correction, on very fine grid (f) CC Mto = 0.7, on very fine grid (1056 x 376) 

(1056 x 376) 


Figure 6. Pressure and heating distributions over a compression corner at Mach 11.3 computed with k — lj 
(2006) turbulence model on refined grids. 
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(e) No compressibility correction, on very fine grid 
(1056 x 376) 

Figure 7. Shear force distributions over a compression corner at Mach 11.3 computed with k — u (2006) 
turbulence model on refined grids. 
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u employs a loosely coupled formulation with lst-order spatial accuracy in the convection terms. Baseline 
LAURA uses strong coupling of all equations so that convective terms are all second-order accurate away 
from any limiting conditions. LAURA was modified so that only the turbulence equations were formulated 
with first-order convection. This change caused a slight decrease in heating levels in the separation zone but 
no change in the extent of separation was observed. The next change to LAURA reset the free stream values 
of k and w to match DPLR; this modification caused no significant change to the pressure or heating. The 
last modification recognizes that baseline LAURA computes the production term V = t 13 (see Eq. 28) 
while baseline DPLR uses V = p, t (fl^ + This modification caused over-prediction of the extent of 

separation and made the result more like DPLR - though it now predicts more separation than DPLR. 

The pressure and heating distributions from three different codes using their standard implementation 
of the Menter-SST model are presented in Fig. 8-b. The LAURA results (black and red) are computed on 
the fine grid. No significant changes are expected if computed on the very fine grid based on the earlier grid 
convergence study in Figs. 6 a versus e, b versus f and in Figs. 7 a versus e and b versus f. The LAURA 
results show a small decrease in the extent of separation for the 5-species air model compared to perfect 
gas. The maximum temperature in the flow is approximately 1,000 K so dissociation is not expected but 
heat capacity includes vibrational modes at these temperatures. The Vulcan results (blue and magenta) are 
computed for a perfect gas on the very fine grid. They show a larger over-prediction in the extent of separation 
for this model as compared to LAURA or DPLR. The Vulcan results were executed with r proi i ii m = 5, more 
restrictive than LAURA with r pro d t um = 20, because unrealistically large values of k and consequently 
are produced in this simulation behind the compression shock if r pro d t nm > 5. Vulcan was rerun without 
any production limiting. This change did reduce the extent of separation though it still was larger than 
indicated by the experimental data and the other simulations. A more comprehensive investigation into the 
causes of these differences has not yet been accomplished. Note that uncertainty metrics for Vulcan runs 
using other models on the same grid are captured in Table 3. DPLR results on the baseline grid for pressure 
are shown in brown. The separation for the baseline SST model in this case is slightly larger than indicated 
by experiment. The LAURA implementation looks good in this case but it does suffer from unrealistically 
large values of k behind the compression shock. 

These differences suggest a new verification test is needed to confirm that model implementation across 
codes is consistent. The method of manufactured solutions could fulfill this need provided that the man- 
ufactured solution was sufficiently complex as to engage all elements of the model (production limiters, 
realizability conditions, etc.). Indeed, the flat plate tests in Appendix B did not reveal significant differences 
in models that were exposed by this wedge problem with separation. A better approach may be to use the 
restart solutions for a complex flow on two successively refined grids, with and without turbulence models 
engaged, as the basis of a “manufactured solution”. The ability to: (1) transfer the restart solutions for 
a complex simulation from the converged solution of one code to another; and (2) follow with a contour 
plot of residuals for every equation prior to any solution update could serve to expose model inconsistencies 
between codes without the need to run the second code to convergence. 



(a) Progression of Models in LAURA 


LAURA -PG P ’ m , Vulcan (limited) 

q, W/m - 



Figure 8. A progression of turbulence model changes and code implementations for the Menter SST model 
applied to Run 54. 
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B. CUBRC - Sharp Cone Flare 


Experimental data for a turbulent boundary layer at Mach 11 on a 6 deg. half-angle sharp cone followed by 
a 36 deg. compression (42 deg. half-angle flare) obtained at the 96 inch CALSPAN shock tunnel 33 are used 
as the second validation check for the suite of turbulence models applied to shock wave turbulent boundary 
layer interactions at compression corners. The simulated configuration is presented in Fig. 9. The simulation 
put a 0.01 inch bluntness on the nosetip. Distances in Fig. 9 are measured from the sphere-cone junction 
point. (Note that gauge locations are recorded 33 as distance from the virtual cone vertex and so appropriate 
transformations are required to compare simulation and experiment.) The baseline grid uses 130 cells along 
the cone and 128 cells across the shock layer. The near wall region prior to separation is resolved with 
y + ft .29, achieved by specifying a cell Reynolds number equal to 0.1 in the initial grid generation process. 
Inflow boundary conditions are: Poo = 1, 807m/s, = 0.032354kg/m 3 , T\ ' x = 67. 4K and a = Odeg.. Surface 

boundary conditions are no-slip with T w = 300K. Supersonic outflow boundary conditions are extrapolated 
from the interior. The test gas, air, is treated as calorically perfect. Molecular transport properties are 
computed from Sutherlands law. 

P 



Figure 9. Cone-flare geometry with shock layer pressures non-dimensionalized by poaV^. 

Pressure and heating distributions for this case using the algebraic turbulence models are presented in 
Fig. 10. The Cebeci-Smith model is in good agreement with extent of separation on the baseline grid (Fig. 
10(a,b)) but significantly over-predicts separation on the next finest grid (Fig. 10(c,d)) which has a factor of 
two increase resolution in all directions. The Baldwin-Lomax model under-predicts the extent of separation 
on both the baseline grid (Fig. 10(e)) and the fine grid (Fig. 10(f)). Good agreement is obtained with 
pressure and heating levels upstream of the interaction and pressure downstream of the interaction using 
the Baldwin-Lomax model. Heating levels are under-predicted by approximately 50% downstream of the 
interaction. The coarse grained metrics from Eqs. 43 - 48 for this figure are captured in Table 4. 

Simulations using the Spalart-Allmaras, Menter SST, and k — to (2006) models are compared to ex- 
periment in Fig. 11. The Spalart-Allmaras model provides best agreement on the baseline grid with the 
experimental data. It has not been tested on the next finer grid. Based on previous grid refinement studies it 
is expected that refinement will cause the extent of separation to further increase. The Menter SST and k — u> 
models without compressibility correction show very little separation on the baseline grid. A compressibility 
correction to the k — io (2006) induces separation on the baseline grid. The coarse grained metrics from Eqs. 
43-48 for this figure are captured in Table 5. 

Refinement studies continue, focusing on the k — u> (2006) models in Fig. 12. Results for M t o = 0.5 
are presented for the coarse and fine grids in Fig. 12 (a,c), respectively. Results for the slightly larger 
argument to the Heaviside function (M t o = 0.7) are presented in Fig. 12 (b,d), respectively. Recall that 
M t o = 0.7 gave the best comparisons to the previous test case (Fig. 6 (b,f)). In this case, the simulation with 
Mf o = 0.7 fails to separate even on the fine grid. (A model without compressibility correction is equivalent 
to engaging the compressibility correction with a value of M t0 that exceeds the turbulence Mach number 
everywhere in the flow. Consequently, the simulation without compressibility corrections fails to separate 
on the fine grid as well.) Best agreement in this case on the fine grid occurs with M t o = 0.5. Pressure 
levels across the recirculation region are well predicted. Both the pressure and heating show an overshoot of 
50% and 80% respectively after reattachment but then recover the measured levels on the flare. (A drop-off 
in agreement beyond x = 109 is not understood at this time but may be a misunderstanding of the gauge 
location definition on the flare.) The coarse-grained metrics from Eqs. 43 - 48 for Fig. 6 are captured in 
Table 6. 
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q, W/cm 2 q, W/cm 2 q, W/cm 




(a) Cebeci-Smith, Eq. 12, (130 x 128) 


(b) Cebeci-Smith, Eq. 10, (130 x 128) 




(c) Cebeci-Smith, Eq. 12, (260 x 256) (d) Cebeci-Smith, Eq. 10, (260 x 256) 




(e) Baldwin-Lomax, (130 x 128) (f) Baldwin- Lomax, (260 x 256) 

Figure 10. Pressure and heating distributions over a 6-degree cone with 36-degree compression at Mach 11 
computed with algebraic turbulence models. 
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,UI/N ‘d UI/N ‘d UI/N ‘d 


Table 4. Coarse-Grained Uncertainty Metrics for Fig. 10 


Model 

E (^Qmax^) 

E(jpmax ) 

E{jlsep) 

E(p S ep ) 

E(c[pre) 

E(j?pre ) 

E(L) 

CS (130 x 128) 

-60 

5 

-130 

-70 

0 

-5 

0 

CSp (130 x 128) 

20 

15 

-10 

10 

60 

-5 

-5 

CS (260 x 256) 

-60 

40 

-90 

-80 

NA 

NA 

NA 

CSp (260 x 256) 

30 

40 

70 

90 

NA 

NA 

-350 

BL (130 x 128) 

-60 

0 

-220 

-130 

-5 

0 

95 

BL (260 x 256) 

-60 

0 

-220 

-130 

-5 

0 

95 


Table 5. Coarse-Grained Uncertainty Metrics for Fig. 11 


Model 

E(qmax) 

E{jPmax ) 

E(q se p') 

E(jp S ep ) 

E^Qpre') 

E(jppre) 

E(L) 

SA (130 x 128) 

-70 

0 

-80 

-200 

20 

-5 

70 

SST (130 x 128) 

160 

0 

-240 

-100 

20 

-5 

100 

k - ui 06 (130 x 128) 

5 

-5 

-240 

-110 

5 

-5 

100 

k — ui 06.25 (130 x 128) 

5 

0 

-240 

-90 

-50 

-10 

40 




(a) Spalart-Allmaras) 


(b) Menter SST 




(c) k — co (2006), no compressibility correc- 
tion 


(d) k — uj (2006), compressibility correction 
M t o = 0.25 


Figure 11. Pressure and heating distributions over a 6-degree cone with 36-degree compression at Mach 11 
computed with 1- and 2-equation models on baseline grid (130 x 128). 
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p, N/m 2 p, N/m 2 


jWS/m ‘b ,uja/M ‘b 


Table 6. Coarse-Grained Uncertainty Metrics for Fig. 12 


Model 

E (Q max ) 

E ( Jpmax ) 

E (Q'sep) 

E {jpsep ) 

E (^pre) 

E ( JPpre ) 

E(L) 

k — uj 06.5 (130 x 128) 

0 

-5 

-240 

-no 

5 

5 

100 

k — uj 06.7 (130 x 128) 

5 

-5 

-240 

-no 

5 

5 

100 

k-uj 06.5 (260 x 256) 

80 

50 

-100 

-5 

5 

5 

-5 

k-uj 06.7 (260 x 256) 

-10 

0 

-240 

-no 

5 

5 

100 




(a) M t0 = 0.50, (130 x 128) 


(b) M t0 = 0.70, (130 x 128) 




(c) M t0 = 0.50, (260 x 256) 


(d) M t0 = 0.70, (260 x 256) 


Figure 12. Pressure and heating distributions over a 6-degree cone with 36-degree compression at Mach 11 
computed with k — uj (2006) model and non-standard compressibility corrections. 
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p, N/m 2 p, N/m 2 


One additional refinement, another factor of two in each direction (520 x 512), is applied to the k — ui 
(2006) model with M t o = 0.5. The results of this simulation on the finest grid are presented in Fig. 13. 
There is a small additional movement of the separation point upstream, seen by comparing Fig. 13 (a) to 
Fig. 12 (c). The pressure and heating overshoots are sharper and larger, now 100% and 80%, respectively. 
Otherwise, the agreement between the fine and very fine grid simulations is good. There is an extremely large 
change in the temperature and velocity profiles going from the baseline (dashed) to the fine (solid) grids in 
Fig. 13 (b). Very small changes are evident in these profiles in going from the fine (solid) to finest (dash-dot) 
grids within the recirculation zone. The near wall region and the inviscid region above the recirculation 
zone show grid convergence on all three grids. The coarse grained metrics from Eqs. 43 - 48 for Fig. 13 are 
captured in Table 7. 




(a) surface pressure and heating, (520 x 512) (b) velocity and temperature profiles at cone- flare junc- 

tion 

Figure 13. Simulation results on finest grid (520 x 512) over a 6-degree cone with 36-degree compression at 
Mach 11 computed with k — u (2006) model and Mto = 0.5. 


Table 7. Coarse-Grained Uncertainty Metrics for Fig. 13 


Model 


E{jpmax ) E(c[sep) E(p S ep') E^fjpre'j ^(.Ppre) £/(Z/) 

k - u> 06.5 (520 x 512) 

100 

80 -90 -20 5 5 -15 


C. Imperial College no. 2 Gun Tunnel, Nitrogen 

Surface pressure 38 and heating 39 measurements for nominal Mach 9.22 flow over a 2D wedge with compression 
angles from 15 to 38 degrees have been used by several authors (see review by Roy and Blottner 12 ) for code 
validation. We capture simulation metrics for several implementations of the k — fl model by Huang and 
Coakley 40 (34 degree ramp) and Coratekin et al. 41 (38 degree ramp) in Table 8. The “res?” notation in 
the table indicates it is not clear that solutions are grid converged in the streamwise direction based on grid 
requirements for the CUBRC wedge tests. The “NA” notation indicates the metric is not available because 
it could not be read from the figures or tables. 

D. Rodi Model Simulations 

A variation of the k — e model by Rodi 42 and tested by Horstman 43 appeared to be one of the better models 
considered in the review by Roy and Blottner. 12 We have not tested the model ourselves but want to cast 
a wide net for the purpose of assessing simulation uncertainty for the SWTBLI at compression corners. 
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Table 8. Coarse-Grained Metrics for Coleman-Stollery Wedge 


Description 


Model 

E ( Qmax ) 

E{jPmax) 

E (Q'sep) 

E{jp S ep) 

E ijipre) 

E (Ppre) 

E(L) 

H&C 


k — u> 88 res? 

80 

-10 

NA 

NA 

NA 

NA 

100 

H&C 

k 

— w 88 cor res? 

0 

2 

NA 

0 

NA 

NA 

-60 

C&vK&B 


k — u> 88 res? 

100 

-20 

90 

-50 

NA 

NA 

40 

C&vK&B 

k 

— to 88 LS res? 

90 

-20 

NA 

NA 

NA 

NA 

40? 

C&vK&B 

k 

— u> 88 RC res? 

-10 

-25 

150 

20 

NA 

NA 

20 


Uncertainty metrics are recorded in Table 9 for simulations of a 2D compression corner (Mach 9.2) , 38, 39 
a cone flare (Mach 9. 2), 33 and the cylinder-flare configuration of Kussoy and Horstman (Mach 7.1) 44 and 
Coleman (Mach 9. 2). 45 Grid convergence tests were reported for these simulations using grids of 80 to 150 
points in the x-direction and 36 to 60 points in the y-direction. These grid resources are significantly less than 
those required in LAURA (see, for example, Fig. 13) for the same simulations to achieve a grid-converged 
solution. 


Table 9. Coarse-Grained Uncertainty Metrics for Rodi Model 


Description 

Mach 

E(qmax ) 

E(Pmax) 

E (Q'sep) 

E^Psep') 

E^Qpre') 

E(jppre ) 

E{L) 

cyl - 20 deg. flare 

7.1 

20 

5 

0 

0 

NA 

NA 

none 

cyl - 35 deg. flare 

7.1 

20 

10 

-40 

20? 

NA 

NA 

0 

2D cc - 34 deg. 

9.2 

35 

-10 

NA 

NA 

NA 

NA 

25 

cyl - 40 deg. flare 

9.2 

35 

-20 

NA 

NA 

NA 

NA 

0 

cone - 36 deg. flare 

11. 

-20 

-10 

NA 

-30 

NA 

NA 

10 


E. CUBRC Compression Corner - Runs 12, 16, and 19 

Free stream conditions for three addition CUBRC Compression corner simulations are provided in Table 10. 
The wedge geometries have a running length of 39.2 in (.99568 m) to the corner. The simulations were run 
with grids of 460 cells in the streamwise direction and 128 cells in the normal direction. Grid metrics are 
believed to be equivalent to the fine grid used in Run 54 simulations. Runs 12 and 16 show no separation at 
the corner and the associated simulations also indicate attached flow. A reduced set of uncertainty metrics 
for these first two cases are presented in Table 11. The experimental data for Run 19 shows separated flow 
but none of the simulations computed any separation. The metrics for this case are captured in Table 12. 


Table 10. Free stream conditions for CUBRC Runs 12, 16, and 19 


Run 

wedge angle, deg. 

Poo, kg/m 3 

Voo, m/s 

Too, K 

Moo 

12 

27 

0.5090 

1390. 

71.1 

8.2 

16 

30 

0.5060 

1375. 

68.9 

8.3 

19 

33 

0.4910 

1379. 

71.7 

8.1 
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Table 11. Coarse-Grained Uncertainty Metrics for Unseparated Flow, CUBRC Runs 12 and 16 



Run 12, 27° 

Run 16, 30° 

Model 

E ((Jrnax') 

E (Pmax ) 

E i^Qmax') 

E ( JPmax ) 

Menter SST 

75 

0 

70 

5 

k — a; 98 

75 

0 

70 

5 

k — a; 06 

50 

5 

45 

0 


Table 12. Coarse-Grained Uncertainty Metrics for CUBRC Run 19 


Model 

E ( Qmax ) 

E (jpmax) 

E (Q'sep) 

E(p S ep ) 

E {j}pre) 

E(Ppre ) 

E(L) 

Menter SST 

30 

-25 

-35 

-75 

15 

0 

100 

k — a; 98 

30 

-25 

-35 

-75 

15 

0 

100 

k — a; 06 

15 

-25 

-35 

-75 

15 

0 

100 


VII. Uncertainty Estimates 

We define a simple uncertainty metric Eqc (CC for compression corner) that tries to capture the diverse 
set of comparison metrics recorded in Tables 2-12. We observe that in current models, if the separation extent 
is well predicted, the post interaction peak values are not well predicted. In like manner, if the post interaction 
peak values are well predicted, the extent of separation is not well predicted. This behavior is captured with 
the following averaged metric, using the most conservative approximation that the experimental uncertainties 
are epistemic - with no knowledge of the probability of where a measurement may exist within the uncertainty 
interval. 

Ecc = Ecc,num + EcC,exp = 7^ {\^(Qmax)\ + \E(p m ax)\ + \E(L)\) + — (\Aq exp \ + |Ap ea; p| + \AL exp \) (49) 

Discretization error is not included here because grid convergence studies indicate discretization error is much 
smaller than the other identified sources of uncertainty. We note that there are fifty-four simulations included 
in the tables involving five experimental data sets from CUBRC, one from Imperial College, and one from 
NASA Ames. There are thirty-six new simulations using LAURA, six new simulations using Vulcan, and 
one new simulation using DPLR. Ten entries in the tables sample simulations that are previously published. 
Table 13 reviews maximum, minimum, average, and median values of Ecc,num for various groupings of 
simulations. 

The first row of Table 13 (All) includes metrics for all models on all configurations. The next eight 
rows of Table 13 include metrics for a specific turbulence model on all configurations where separation 
was observed experimentally. Of these models, the sixth row (fc - w 06 - cmp) includes compressibility 
corrections with various values of M t o; other compressibility corrections and other techniques for turning 
off their influence near the wall have not been tested. The tenth row (unseparated) gathers metrics for 
simulations of configurations where separation was not observed experimentally. The simulations tend to 
do a better job for attached flows and so these cases were isolated so as not to skew results for the more 
challenging separated flow cases. The last row (fine, very fine) only considers simulations that were done on 
the two finest grids for which discretization error is demonstrably smaller than other sources of error. 

In computing a representative metric of uncertainty for production ready models as discussed in the 
introduction we discount some entries in Table 13. Results from the algebraic models (row 7) are not included 
because these models failed to converge for several separated flow conditions. Results from earlier simulations 
in rows 8 and 9 are not included because of concerns that the grid densities used to compute these cases were 
considerably less than grid densities required here to demonstrate convergence. Continued grid refinement 
to the fine and very fine levels used herein can actually degrade agreement with experimental data for some 
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Table 13. 


norms 


Set Description 

Number in Set 

EcC,num,min 

EcC,num,max 

E(J C,num, avg 

EcC,num,med 

Notes 

All 

53 

8.3 

140.0 

43.5 

40.0 

Omit 4.3 

SA 

2 

20.0 

46.7 

33.3 

33.3 


SST 

7 

28.3 

106.7 

65.2 

51.7 

Omit unsep. 

k — u> 98 

4 

40.0 

55.0 

49.6 

51.7 

Omit unsep. 

k — u> 06 

7 

33.3 

48.3 

42.1 

43.3 

Omit unsep. 

k — 10 06 - cmp 

12 

15.0 

86.7 

49.9 

41.7 

Sarkar/Zeman 

Algebraic 

5 

13.3 

140.0 

55.7 

51.7 

Omit 4.3 

k — cj 88 

5 

18.3 

63.3 

41.1 

50.0 

Table 8 

Rodi 

4 

10.0 

23.3 

16.3 

15.8 

Omit unsep. 

unseparated 

7 

8.3 

26.7 

20.5 

25.0 

various models 

fine, very fine 

15 

33.3 

106.7 

60.0 

55.0 

various models 


models because the post shock peaks sharpen and rise. (See, for example Fig. 4-e versus Fig. 6-d and Fig. 
12 (a versus c).) Consequently, the metrics in the last row of Table 13, specifically Ecc,num,med = 55, are 
thought to be the most representative assessment of the numerical component of uncertainty for simulations 
of SWTBLI at compression corners sufficiently steep to cause separation. In like manner, Ecc,num,med = 25 
from row 10 is considered a reasonable assessment of the numerical component of uncertainty for simulations 
of SWTBLI at compression corners where flow remains attached. Using experimental uncertainty estimates 
from Sec. VI - B the total uncertainty for separated flows (or conditions near incipient separation) is 

Ecc = E cc ,num + E cc ,ex P = 55 + (9.1 + 5.8 + 12.6)/3 = 64 (50) 

The corresponding estimate for attached flows, away from an incipient separation domain, is 

Ecc = Ecc,num + Ecc, exp = 25 + (9.1 + 5.8 + 0)/3 = 30 (51) 

Clearly, this study has not evaluated every permutation of turbulence model and associated corrections 
to provide the best possible value of Ecc ■ Before adding new simulations it is critical that a well defined 
verification process be established so that consistency of models across codes can be certified. For now, 
we have defined a process involving a collection of test cases to provide an estimate for the state-of-the-art 
using standard models. We expect this metric to serve as a reasonable goal to evaluate improvements to the 
state-of-the-art for compression corner problems. As simulations improve, a more comprehensive, perhaps 
fine-grained metric, will be required. For example, the current definition of Ecc,num could equal zero and 
still have significant errors in the plateau region predictions for heating and pressure. This simple metric and 
target uncertainty serve the original purpose of providing the Fundamental Aeronautics Program with useful 
parameters and processes for evaluating improvements to the state-of-the-art in hypersonic simulation. 

VIII. Simulations of Mission Relevant Test Problem 

The mission relevant test problems were defined previously in Sec. II. New simulations are executed 
to assess additional uncertainties due to high temperature effects. Higher temperatures are progressively 
encountered as Mach number is increased and as compression angle is increased. The Spalart-Allmaras (SA) 
model and Menter-SST model have been applied to the test problems with perfect gas and 5-species models 
for air specifying equilibrium, finite-rate (chemical nonequilibrium), and frozen flow. Representative results 
follow. 

The Mach 7 simulations indicate peak temperature in the boundary layer approaching the interaction 
is approximately 780 K. The post-compression peak temperature is approximately 1,260 K (nonequilibrium 
model) or 1,380 K (perfect-gas model). There is no significant chemical dissociation at these conditions. 
Consequently, all of the 5-species models are expected to yield equivalent results because they use identical 
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thermodynamic and transport property models. Figure 14 for Mach 7 confirms this expectation. It shows 
that equilibrium and non-equilibrium results are nearly identical on the scale of the figure for shallow (a) 
and steep (b) compression angles using two different turbulence models. The perfect-gas model shows a 
14% high heating level compared to the high temperature gas models for the shallow compression with the 
SA model. The perfect gas pressure distribution is in excellent agreement with the high temperature gas 
model pressure distributions. At Mach 7 it is evident that the perfect-gas models introduce some error in 
heating distribution due to use of constant heat capacity and use of Sutherlands law. The augmentation of 
uncertainty at Mach 7 due to high temperature gas effects for the mission relevant test problem is therefore 
negligible relative to the baseline model form uncertainty defined previously assuming an appropriate gas 
model is used. 


Equil. 




(a) Mach 7, 5.5° compression, SA model, (500 x 64) (b) Mach 7, 30° compression, SST model, (500 x 64) 


Figure 14. Effects of gas model on simulation of surface pressure (black lines) and heating (red lines) for 
turbulent flow over test configuration at Mach 7. 


Figure 15 shows shows surface pressure and heating levels across the interaction region of the 5.5° com- 
pression test case using the Menter SST model at Mach 14. The pressure levels approaching the interaction 
are essentially independent of the 5-species gas model used in the simulation. The non-equilibrium simu- 
lation shows a 15% higher heating than the frozen or equilibrium chemical models. The difference here is 
caused by the use of a fully catalytic boundary condition in all simulations though only the non-equilibrium 
simulation has a significant catalytic (diffusion) component. In the frozen flow model there is no dissociation 
so there is no catalytic component to heating. In the equilibrium flow model the temperature in the near 
wall region is too low to enable any significant level of atoms; so catalytic heating is again suppressed. The 
non-equilibrium simulations indicate peak temperature in the boundary layer approaching the interaction is 
approximately 2,270 K and the oxygen mass fraction is approximately 0.054, mostly coming from a blunted 
leading edge upstream as will be discussed subsequently. In this case the atoms diffuse to the surface and a 
catalytic component to heating is realized. 

There is no clearly defined pressure or heating maximum in the post interaction domain. This case is 
unique in that it was run using an earlier version of the mission relevant problem configuration with a 1 
cm leading edge bluntness. Early tests on this configuration with laminar flow (not shown here) indicated 
entropy effects from streamlines passing through the curved, leading edge bow shock persist to the interaction 
and cause a more gradual rise in pressure and heating after the compression. This issue is more pronounced 
for small compressions and large incoming boundary layers. It is evident to a lesser extent in the earlier 
Mach 7 result of Fig. 14 (a) even though that simulation had a sharp leading edge. 

Progressing now to the 30° compression at Mach 14 the post-compression temperature in the shear layer 
behind the bow shock - compression shock interaction is approximately 3,200 K and significant oxygen 
dissociation follows in the free shear layer from that crossing. Closer to the surface, the temperature behind 
the compression shock near the wall approaches 3,600 K with atomic oxygen mass fraction approaching 
0.0016. Note that this atomic oxygen mass fraction is considerably less than noted in the previous case. 
In the previous example with the blunted leading edge significant levels of atomic oxygen produced at the 
blunted leading edge convect to the compression corner. In the current case the reservoir of atoms that 
drives catalytic heating has not been augmented from sources at the sharp leading edge. Because there 
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is relatively little dissociation near the wall across the compression corner one might expect insignificant 
differences between 5-species model simulations using chemically frozen, equilibrium, and non-equilibrium 
gas models. Indeed, the catalytic heating component is less than 0.05% of the heating rate just upstream of 
the interaction and distributions of pressure and shear for all 5-species models differ by less than 5% (Fig. 
16). Only the perfect-gas model, inappropriate for these temperatures, shows significant differences (7% for 
pressure, 20% for heating). Uncertainty metrics should not be modified based on use of an inappropriate 
gas model at these mission relevant conditions. 




(a) pressure (512 x 128) (b) heating (512 x 128) 

Figure 15. Effects of gas model on simulation of surface pressure (black lines) and heating (red lines) for 
turbulent flow over test configuration with 5.5° compression at Mach 14 using Menter SST model. 




(a) pressure (500 x 64) 


(b) heating (500 x 64) 


Figure 16. Effects of gas model on simulation of surface pressure (black lines) and heating (red lines) for 
turbulent flow over test configuration with 30° compression at Mach 14. 


Simulations of the 5.5° and 30° compression corner configurations did not show any separation at Mach 
7 or at Mach 14. Simulations of a 35° compression corner problem using the Menter SST turbulence model 
with non-equilibrium gas chemistry have been executed on several different grids in the interaction region 
with focus on resolving an incipient separation. Two of these simulations are considered in Fig. 17. Grid 4 
uses the standard grid distribution function and shock alignment algorithms as described earlier for CUBRC 
Run 54. Grid 2 adjusts the shock alignment algorithm so that the distance between the wall and the inflow 
boundary downstream of the corner is never allowed to be less than the distance computed by the alignment 
algorithm at the corner. Grid 4 enables a tighter resolution of the compression shock at the expense of 
significant grid skewness. Grid 2 has much less grid skewness but the compression shock cuts the grid along 
diagonals. An overlay of the results using Grid 2 and Grid 4 show good agreement of shock shape and 
pressure contours. The extent of separation is still small, extending only 2 cells (sa 3 mm) upstream from 
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the corner on the finest grid. 

Corresponding pressure and heating distributions (Fig. 18) show small differences between the peak 
values (less than 5% on pressure and less than 8% on heating) for the same models on two fine grids. The 
perfect gas results again show more significant differences indicating they are not appropriate for this flow 
condition. 




(a) Grid topologies colored by log(pressure) 


(b) Overlay of pressure contours 


Figure 17. Grid topologies used in simulation of test configuration with 35° compression at Mach 14. 




(a) pressure (b) heating 

Figure 18. Effects of gas model and grid on simulation of surface pressure (black lines) and heating (red lines) 
for turbulent flow over test configuration with 35° compression at Mach 14 using Menter SST model. 


The following observations are offered to summarize the impact of chemistry to the uncertainty metrics 
for mission relevant problems with a sharp leading edge. 

• At Mach 7 there is no significant production of atoms and any 5-species air model gives equivalent 
results. 

• At Mach 14 dissociation of atomic oxygen comes in to play, especially for the steeper compressions. 
The greatest difference in any uncertainty metric as a function of 5-species model used is less than 5%. 
This value is considered an upper bound to gas chemistry component of uncertainty. 

• Perfect gas models introduce some error at Mach 7 and Mach 14 and are generally not appropriate for 
mission relevant problem simulations. A perfect gas model shows differences (errors) associated with 
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its assumption of constant specific heats at temperatures above 1,000 K encountered at Mach 7 and 
above. 

• Nose bluntness may introduce additional levels of uncertainty associated with chemistry due to the 
convection of atoms from the nose to the interaction region. The catalytic efficiency of the wall and the 
distance from the nose to the compression corner become additional factors that impact the uncertainty 
metrics. 


IX. Concluding Remarks 

An uncertainty metric is defined for simulation of a shock wave turbulent boundary layer interaction 
(SWTBLI) at two-dimensional or axi-symmetric compression corners using production-ready turbulence 
models. The metric is intended for use by NASA’s Fundamental Aero Program to provide a baseline circa 
2006 for judging improvements to hypersonic simulation as a result of investments in modeling and experi- 
ments. 

A mission relevant test problem is defined including free stream conditions for which the simulation 
uncertainty will be characterized. The governing equations and turbulence models are fully documented 
herein to remove any ambiguity of how models are formulated. A set of coarse grained uncertainty metrics 
are introduced to characterize the essential elements of a SWTBLI within a simulation. These metrics 
are also defined to easily capture their values from figures in the literature. Relevant experimental data 
on compression corners are reviewed and new simulations are executed with a variety of models and codes 
representative of the state of the art in 2006. The metrics are calculated for both existing and new simulations 
and a median value equal to 55% is computed for SWTBLI with separation on the basis of new, grid converged 
simulations. It is observed that grid converged solutions may enable sharper peaks in post shock pressure 
and heating as well as a larger extent of separation that may degrade comparisons with experimental data. 
Finally, new simulations are executed on the mission relevant problem to include gas chemistry perturbations 
in order to assess an overall uncertainty metric. Two appendices are also provided to present verification 
and validation data on a simpler flat plate problem. 

Two critical needs are identified if further reductions in the uncertainty metrics are required. The first 
need, verification of models between codes, is an exceptionally challenging problem for separated SWTBLI 
simulations. A verification exercise on flat plates between LAURA and DPLR failed to identify an implemen- 
tation difference in the models which was critical in the SWTBLI simulation. In like manner, comparisons 
of simulations using Vulcan with LAURA and DPLR show differences which are thought to follow from 
some subtle difference in model implementation but as yet are not explained. It is suggested that the restart 
solutions for a complex flow on two successively refined grids from a single code could be used as the basis of 
“manufactured solutions” for other codes. Plotting the residuals from the second code after a single iteration 
starting from the “manufactured solution” could expose model inconsistencies between codes without the 
need to run the second code to convergence. The second need, refinement of experimental error estimates, 
requires both additional experiments covering relevant energies and Reynolds numbers as well as experiments 
devoted to better define uncertainties in calibration, free stream conditions, flow non-uniformities (temporal 
and spatial), and boundary conditions. 
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Appendix A: Analytic Approximations 


Van Driest Transform I 

Huang and Coleman 46,47 investigate the generality of the Van Driest transformation I 48 on DNS simulations 
of compressible, wall bounded flows. They found that the transformation tended to collapse the computed 
velocity profiles to match the incompressible logarithmic law for cold wall channel cases at Mach numbers 
of 1.5 and 3. In contrast to DNS simulations the velocity profiles computed with RANS models are derived 
from these fundamental relations. Consequently, transformed velocity profiles provide no new fundamental 
insight regarding the physics of turbulent boundary layers. Nevertheless, the transformed velocity profiles 
for the hypersonic simulations are computed to see if the collapsing trends observed by Huang and Coleman 
are verified in current models. 

Given a computed or measured velocity profile U + (y + ) the Van Driest transformed profile Uy D (y + ) as 
presented by Huang and Coleman is given by 


where 


Uvd ~ 


Rvd 


sin 


Rvd{U + + Hvd ) 
Dvd 


— sin 


Rvd Hvd 
Dvd 


Rvd = M r yJ (7 — l)Pr t /2 
Hvd = B q /[{ 7 -l)M 2 ] 
Dvd = \/l + (Rvd Hvd) 2 - 


(52) 


(53) 

(54) 

(55) 


The transformed profiles in the RANS simulations are based on computed values of r w and q w . Note that 
B q and therefore Hvd are negative quantities for cold walls (energy flux into wall, opposite increasing y + 
direction) . 


Van Driest Transform II 

The Karman-Schoenherr equation, 49-51 relating skin friction to momentum thickness Reynolds number for 
an incompressible flow, is given by 

C f = [I7.08(log 10 Reg) 2 + 25.11(log 10 Reg) + 6.012] _1 (56) 

Given a numerical simulation of a compressible flow with computed distribution of Reg a mapping to a 
corresponding incompressible value Reg to be used in Eq. 56 may be implemented with 

Reg = FgReg. (57) 


The incompressible skin friction coefficient can then be computed in Eq. 56 and a transform back to the 
compressible value is implemented with 

C f = C f /F c . (58) 

The transformation factors Fg and F c were derived by van Driest 52 and are expressed as presented by Hopkins 
and Inouye 51 by 

Fg — Me/ Mu; (59) 


and 


where 


and 


olvd = 


F c = C V d [sin 1 {a VD ) + sin 1 (@vd)\ 
2Ay D — B V d a B V d 


-1 


\f^Ay D + By D 


> 0VD — 


\f^Ay D + By D 


a vd = 


CypT e 

T„, 


Byo = 


T e + CypT e - T w 
T,„ 


Cyp = 0 .qA-±M 2 


(60) 

(61) 


(62) 


Comparisons of computed skin friction using turbulence models to the van Driest II transforms are imple- 
mented as a supplemental verification check. 
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Appendix B: Verification and Validation on Flat Plate at Mach 11 

Experimental data for a turbulent boundary layer at Mach 11.1 on a flat plate obtained at CUBRC 37 
are used as a preliminary verification and validation of the suite of turbulence models. Implementation of 
turbulence models is fraught with opportunities to err in the programming of the model or the coupling of 
the model with the conservation equations. A verification check is made by comparing to incompressible 
profiles from the Van Driest transformations discussed in Appendix A. 

X 52” X 

^ 

• Flat Plate only 

• 0° angle of attack 

Figure 19. Flat plate geometry for CUBRC Run 7. 

The geometry is presented in Fig. 19. The baseline grid uses 300 cells along the plate and 96 cells 
normal to the plate. The near wall region is resolved with y + « .035. Inflow boundary conditions are: 
Voo = 1, 780m/s, poo = 0.09483kg/m 3 , T; ' x = 64K and a = Odeg.. Surface boundary conditions are no-slip 
with T w = 300K. Supersonic outflow boundary conditions are extrapolated from the interior. The test gas, 
air, is treated as calorically perfect. Molecular transport properties are computed from Sutherlands law. 

Figure 20 compares experimental data for shear and heating along the surface to simulations using various 
turbulence models. The solid line result in each figure shows the k — ui (2006) model distribution - - providing 
a visual reference using the most recent model of the group. The dashed line provides results for the other 
models. A laminar result is presented in Fig. 20 (a) as a dashed line. The laminar heating result agrees with 
experimental data for heating prior to transition. The turbulent pressure (black line) exceeds the laminar 
value due to a larger displacement thickness. 

The solid line, k — to (2006) without compressibility correction, agrees well with experimental data for 
shear but over-predicts the turbulent heating rate. (See Sec. V for a discussion of associated experimental 
uncertainties.) The algebraic models (Fig. 20 (b,c)) show shear levels slightly less than the reference k — to 
(2006) model but still in agreement with the spread in experimental data. The algebraic model heating levels 
are in better agreement with the experimental data though they still tend to over-predict the heating levels. 
The 1-eq. Spalart Allmaras model (Fig. 20 (d)) is in good agreement with the k — to (2006) model without 
compressibility correction. There is no pseudo transition point in either the algebraic or 1-eq. models - an 
increased level of heating and shear exceeding the laminar value persists to the leading edge. The Menter 
SST model (Fig. 20 (e)) and the older k — u (1998) model (Fig. 20 (f)) are bounded by the predictions 
from the algebraic models and the k — to (2006) model. All models are in reasonably good agreement with 
experimental data for shear and all models exceed the experimental data for heating rate. 

The Van Driest II transform (see Appendix A, section B) is used as a verification check of the shear 
distribution on a flat plate for all models in Fig. 21. The solid red line in these figures is the computed 
momentum thickness Reynolds number from the simulation. The dashed red line is the transformed incom- 
pressible momentum thickness Reynolds number. The solid blue line is the shear distribution calculated 
from the simulation. The dashed blue line is the Van Driest II shear distribution which is obtained by 
assuming the Karman-Schoenherr (Eq. 56) can be used to plot the incompressible shear as a function of 
the transformed momentum thickness Reynolds number and that Eq. 58 can be used to transform the 
incompressible shear back to a compressible distribution. The shear distribution computed with the k — u> 
(2006) model without compressibility correction (Fig. 21 (a)) is in excellent agreement with the Van Driest 
II transform distribution. The algebraic models (Fig. 21 (b,c)) tend to produce a slightly lower value Reg 
and a correspondingly higher value of transformed shear. Nevertheless, both values of shear are in good 
agreement with experimental data. The 1-eq. Spalart-Allmaras model (Fig. 21 (d)) and the 2-eq Menter 
SST models (Fig. 21 (e)) slightly under-predict the Van Driest II levels but remain in better agreement 
than either of the algebraic models. The older k — to (1998) (Fig. 21 (f)) appears equivalent to the algebraic 
models in terms of the level of agreement with the transformed shear. 
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The Van Driest I transform (see Appendix A, section A) serves as a second verification check by comparing 
the transformed, incompressible velocity profile to the classical incompressible result with u + = y + in the near 
wall (inner) region and u + = In (y + )/ n + C for y + > 10 to the boundary layer edge (outer). The transformed 
velocity profiles for each of the turbulence models is presented in Fig. 22 where the compressible velocity 
profile, non-dimensionalized by u T is given by the solid black line, the inner and outer reference profiles 
with k = 0.41 and C = 5.0 are given by the dashed red and blue lines, respectively, and the Van Driest 
transformed velocity is given by the dash-dot red line. All of the models follow the inner law variation for 
y + < 3. The k — u (2006) model without compressibility correction (Fig. 22 (a)) and the Menter SST model 
(Fig. 22 (e)) have transformed velocity profiles that are in excellent agreement with the classical limit in the 
outer region ( y + > 10). The other models remain in fair agreement with the outer limit — small changes 
in the value of n or C would produce better agreement in these cases. All of the models are in acceptable 
agreement with the classical limit in this verification check. 

The effects of the compressibility correction (see Eq. 29) in the 2-eq. models are explored in Figs. 23- 
25. The compressibility correction is not intended for wall bounded flows. Rather, it is applied to problems 
involving mixing in free shear layers. The thought here is that compression corners will often induce separated 
flow and the free shear layer in these cases may not be adequately simulated without some compressibility 
correction. The concern is that the compressibility correction may significantly distort the results achieved 
for the fully attached case on the flat plate. Consequently, these corrections are tested here. 

Figure 23 follows the same format as Fig. 20. It compares experimental data for shear and heating along 
the surface to simulations using various turbulence models with a compressibility correction. The solid line 
result in each figure again shows the k — u (2006) model distribution without compressibility correction. The 
dashed line provides results for the other models now including a compressibility correction. The baseline 
compressibility correction with M t o = 0.25 applied to the k — u (2006) model is shown with the dashed lines 
in Fig. 23 (a). The heating in this case now slightly under-predicts the data and the shear, which had been 
in good agreement with experiment, now significantly under-predicts the data. As the critical turbulent 
Mach number is increased (from M t0 = 0.3 in Fig. 23 (b) to M t0 = 0.5 in Fig. 23 (d)) the heating and 
shear distribution approach the compressibility correction free model. A value of 0.3 < M t o < 0.4 slightly 
compromises the agreement of the baseline model with shear data but improves the agreement with heating 
data. The compressibility correction behavior for the Menter SST model (Fig. 23 (e)) and the k — u (1998) 
model (Fig. 23 (f)) behave similarly to the k — u (2006) model. 

The compressibility correction with M t0 < 0.3 abrogates agreement noted earlier with the Van Driest 
II transformed shear distribution in Fig. 24 and with the Van Driest I profiles in Fig. 25. Only the cases 
with M t o = 0.4 (subfigure (c)) and M t o = 0.5 (subfigure (d)) show agreement with the incompressible 
distributions. Disagreement here is inconsistent with existing DNS simulations at lower edge Mach numbers. 
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(a) Laminar (b) Algebraic: Baldwin-Lomax 




X 


(c) Algebraic: Cebeci-Smith. 


(d) 1-Eq: Spalart-Almaras (Catrix-Aupoix Com- 

pressible Model) 




(e) 2-Eq: Menter SST 


(f) 2-Eq: k-u (1998) 


Figure 20. Pressure, shear, and heating distributions over a flat plate at Mach 11.1 computed with the k — lj 
(2006) model compared to experimental data and simulations with other models. 
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(a) 2-Eq: k - u (2006) 



(c) Algebraic: Cebeci-Smith. 



(b) Algebraic: Baldwin-Lomax 



(d) 1-Eq: Spalart-Almaras (Catrix-Aupoix Compress- 
ible Model) 




(e) 2-Eq: Menter SST 


(f) 2-Eq: k - uj (1998) 


Figure 21. Shear distribution computed from turbulence models compared to shear computed from Van Driest 
transformation. The computed momentum thickness Reynolds number distribution from each model and the 
transformed, incompressible momentum thickness Reynolds number used in the Van Driest transform show 
consistent magnitudes for all models. 
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(b) Algebraic: Baldwin-Lomax 



(c) Algebraic: Cebeci-Smith. (d) 1-Eq: Spalart-Almaras (Catrix-Aupoix Compress- 

ible Model) 




Figure 22. The U+ profile for various models (black line) at x = lm transformed using Van Driest theory to 
an incompressible profile (brown dashed dot line) with comparisons to the incompressible limits for y+ < 1 
(red dashed line) and y + > 10 (blue dashed line). 
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(a) k-w (2006), M t0 = 0.25 (b) k-ui (2006), M t0 = 0.30 




(c) k-u (2006), Mto = 0.40 (d) k-w (2006), Mto = 0.50 




(e) Menter SST, Mto = 0.25 (f) k — u (1998), Mto = 0.25 


Figure 23. Effects of Wilcox compressibility corrections on pressure, shear, and heating distributions over a 
flat plate at Mach 11.1 compared to experimental data. 
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(p/10) N/m 2 T x , N/m 2 (p/10)N/m 2 t x , N/m 2 (p/10)N/m 2 t x , N/m : 




(a) k-cj (2006), Mto = 0.25 (b) k - uj (2006), M t0 = 0.30 




(c) k-Lj (2006), M t0 = 0.40 


(d) k-uj (2006), Mto = 0.50 




(e) Menter SST, Mto — 0.25 


(f) k-w (1998), Mto = 0.25 


Figure 24. Shear distribution computed from turbulence models using Wilcox compressibility corrections 
compared to shear computed from Van Driest transformation. The computed momentum thickness Reynolds 
number distribution from each model and the transformed, incompressible momentum thickness Reynolds 
number used in the Van Driest transform show sensitivity to choice of turbulent Mach number cutoff parameter 

M t o- 
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(a) k-w (2006), M t0 = 0.25 


(b) k-u (2006), M t0 = 0.30 



(c) k-uj (2006), M t0 = 0.40 


(d) k-w (2006), M t0 = 0.50 



(e) Menter SST, Mto = 0.25 


(f) k-u (1998), M t0 = 0.25 


Figure 25. Effect of Wilcox compressibility correction on the 17+ profile for various models (black line) at x = lm 
transformed using Van Driest theory to an incompressible profile (brown dashed dot line) with comparisons 
to the incompressible limits for y+ < 1 (red dashed line) and > 10 (blue dashed line). Agreement with the 
incompressible limit for y + > 10 is degraded as the turbulent Mach number cutoff Mto is decreased. 
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Sensitivity of the simulated results to Prt is investigated in Fig. 26. The dashed lines in each subfigure 
represent the new results with Prt = 1.0 (a) and Pr t = 0.5 (b). A smaller value of Pr t corresponds to a 
larger value of turbulent conductivity. In the case of the flat plate, the temperature peak in the boundary 
layer is diffused and the surface heating is reduced when Pr t is decreased. (See Fig. 27 (a).) In contrast, 
note that a lower value of the critical turbulent Mach number in the compressibility correction (Fig. 27 
(b)) tends to move the temperature peak further away from the wall and lower the heating rate. In general, 
agreement with experimental data for heating is improved but agreement with shear is worse with Pr t =0.5 
for this case. 




(a) Prt = 1.0 (b) Prt = 0.5 

Figure 26. Effects of turbulent Prandtl number Prt on pressure, shear, and heating distributions over a flat 
plate at Mach 11.1 in the k — w (2006) model. 




(a) Prt effects with no compressibility correction 


(b) Compressibility correction effects with Prt = 0.9 


Figure 27. Temperature profile at x = lm as a function of turbulent Prandtl number and compressibility 
correction for Mach 11.1 flow over a flat plate in the k — uj (2006) model. 
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The next two sets of figures on the flat plate simulations address issues of turbulent kinetic energy coupling 
(Fig. 28) and grid convergence (Fig. 29). In the present formulation turbulent kinetic energy, fc, is included 
as part of the total energy, (i.e. E = e + ( u 2 + v 2 + w 2 )/ 2 + k). Consequently, the convective and diffusive 
transport of k are included in the total energy conservation, Eq. 3. If one assumes that the turbulent kinetic 
energy content of the flow is significantly smaller than the kinetic and internal energy content everywhere 
then an uncoupled energy formulation may be tested with all terms involving k omitted from Eq. 3 and 
E = e + ( u 2 + v 2 + w 2 )/2. Plots of surface shear and heating in Fig. 28 (a) and the temperature profile 
at x = lm in Fig. 28 (b) with default coupled results in solid lines and uncoupled results in dashed lines 
show negligible effect of decoupling in this application. The same set of distributions are used to show grid 
convergence in Fig. 29. The baseline grid of (300 x 96) cells was decreased by a factor of 2 in both directions 
to (150 x 48) cells. There is almost a perfect overlay of results on the two grid systems. 




(a) Surface distributions 


(b) Temperature profile at x = lm 


Figure 28. Effects of turbulent kinetic energy coupling on simulations of Mach 11.1 flow over a flat plate in 
the k — uj (2006) model. 


Finally, we compare results from implementation of models in LAURA and DPLR (V4.01.1) in Fig. 30 
using laminar, algebraic Baldwin-Lomax, Spalart Almaras, and Menter-SST models. Original implementa- 
tions of the models are shown on the left (a) and more consistent implementations of the models are shown 
on the right (b). Simulations using laminar and algebraic turbulence models are in excellent agreement 
between the two codes. Small differences are noted between the 1- and 2-equation model results as originally 
implemented in the two codes. The difference in the 1-equation model came from an extra source term as- 
sociated with different interpretations of the Catrix-Aupoix compressibility changes to the Spalart- Almaras 
model. The original presentation of the model was not written in strong conservation form. LAURA added 
a source term to enable the equation to be written in strong conservation form. DPLR coding did not alter 
the source terms in this way. When this extra source term was removed (resulting in Eq.20 as the new 
baseline) excellent agreement between implementations was achieved. The difference in the Menter-SST 
simulations was traced to use of lst-order accurate implementation of convective flux for the two turbu- 
lence equations in DPLR (loosely coupled) and 2nd-order accurate implementation of convective flux for all 
equations in LAURA (strongly coupled). When the LAURA implementation was modified to use lst-order 
flux the agreement is much improved. Note that differences in the implementation of the production term 
(discussed in association with Fig. 8) produced changes that were small compared to the difference in order 
of accuracy for this flat plate case. The laminar simulation is in good agreement with experimental data 
up to the transition location. After transition the algebraic turbulence model provides best agreement with 
experimental data for the flat plate. 
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wo/M 




(a) Surface distributions (b) Temperature profile at x = lm 


Figure 29. Effects of grid resolution on simulations of Mach 11.1 flow over a flat plate in the k — co (2006) 
model. 




(a) Original models in LAURA and DPLR 


(b) Fully consistent models in LAURA and DPLR 


Figure 30. Comparison of turbulence model simulations as implemented in LAURA and DPLR. 
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